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Abstract The aligned pulsars whose rotation axes and magnetic dipole axes are 
parallel should be positively charged. The total charge of pulsars is calculated after 
considering the electromagnetic field in and out the star under a specific condition. 
The statistical relation between the pulsar's rotation energy loss rate (or the period 
derivative) and the period may hint that the millisecond radio pulsars with small 
periods could be low-mass bare strange stars. 
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1 INTRODUCTION 

The study of pulsar magnetosphere is essential for us to understand various radiative processes, 
and thus observed emission in different bands. The charge-separated plasma (Goldreich and 
Julian 1969) in the magnetosphere of various gap models (Ruderman & Sutherland 1975; Arons 
& Scharlemenn 1979; Cheng, Ho, & Ruderman 1986; Qiao et al. 2004) has been investigated. 

After investigating the electromagnetic field of an aligned pulsar whose rotation axis is 
parallel to the magnetic axis, Xu, Cui & Qiao (2006) suggested that these pulsars should be 
positively monopolc-charged. We would like to study further charged pulsars in this paper. 
Supposing the pulsar is magnetized homogenously, we calculate the total charges Q. Provided 
that this charge is proportional to the radiation of pulsars, Q should be proportional to E, 
the spin down energy loss rate. The statistical relation between the pulsar's E (or the period 
derivative P) and the period may hint that the millisecond radio pulsars with small periods 
could be low-mass bare strange stars. 

This paper is arranged as follows. The total charge of pulsars is estimated in §2. In §3, we 
show possible evidence for low-mass millisecond pulsars with small P by observational data. 
Conclusions are presented in §4. 

2 PULSAR'S CHARGE ESTIMATED 

If the pulsar is assumed to be magnetized homogenously, its magnetic field could be described 

as 

Bin = = — (r < ro), (1) 
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and 



^^^^^^ 3rir.rn)^m ^ ^^^^ 



where M is the magnetized intensity, m is the magnetic dipole moment, ro is the radius 
of pulsar, r is the distance from one arbitrary point in space to the stellar core, r is the 
corresponding unit vector. Here we define l-Binj = Bq. 

If all the particles in and out the pulsar rotate synchronously around the magnetic axis, the 
electric field, E, could satisfy 

E+ xB^O, (3) 

c 

without considering the magnetic field induced by synchronously rotating particles, where 
is the angular velocity of star rotating around the dipole rotation axis, which relates with 
rotating period of star P by P = 2Tr/n. If k is the unit vector of angular momentum, — Qk 
and m ~ mk for aligned pulsars. We estimate the charge under this specific condition. 
From Eqs.lQJ to lO, the electric fields in and out a pulsar are, respectively, 

^^in = X (fc X f) (r<ro), (4) 

c 

E out = ^^[Hk-f)f~k]x{kxf) (r>ro). (5) 

The corresponding charges then could be 

_ 4^rg _ W-Ein 4^rg _ 2rgr>Po 

3 " 4^ 3 " 3c ' ^' 

and Qout = / Pout • dV = 0, where Pin and Pout are the charge densities in and out the star. 
From Eqs.l^J and ||SJ), the charge intensity and the charges on the stellar surface is 

Cs — -; — (,-C'out ~ ^in)r=ra ' f — 7, j {' ) 

477 Sttc 

a^-ds^-^ , (8) 

r—ro ^ 

where 9 is the polar angular between arbitrary vector r and k. 
From the calculation above, we can obtain the total charges 

gtotai ^ Qin + Qout + Qs^ ^^fj^" ^ 2.3 X IQioIlgH (coulombs), (9) 

where rg = 7-0/(10^ cm), B12 — B^/ilQ^^ G). From this equation, it is shown that the aligned 
pulsars should be positively charged. 

It is evident that this total charge in Eq. lO is different from that obtained in Eq. (17) of Xu 
et al. (2006, Q ~ 10~^r|i?i2/P^). This difference could result in a shift of the boundary (i.e., 
the field lines being at the same electric potential as the interstellar medium) which separating 
regions I and II in Fig. 1 of Xu et al. (2006): the boundary should be closer (i.e., above the 
critical field lines) to the magnetic axis if the real charge is greater than Q. Unfortunately, we are 
faraway from obtaining an exact value of charge since a global solution to pulsar magnetosphere 
is still impossible. 
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Fig. 1 Left upper panel: Period derivatives P as a function of period P for 1349 

radio pulsars. The dots represent the normal radio pulsars with B > Bq and the 
stars are the millisecond radio pulsars with B < Bq if the dividing line of them is 
the magnetic field Bq = 5 x 10^° G. Left lower panel: The period derivative P as a 
function of period P for millisecond radio pulsars. The empty circles (Group /') and 
squares {Group II) are two groups divided by Kmeans methods. The regions between 
two solid lines(for Group I') and dash lines {Group W) represent limit ranges of radii 
for these two groups, respectively. Right upper panel: Spin down energy loss rate E 
as a function of period P for millisecond radio pulsars. The empty circles {Group T) 
and squares {Group IT) are two groups divided by Kmeans methods. The solid and 
dash lines represent the minimums of radius for these two groups. Right lower panel: 
Two-peak period P distribution of 87 millisecond radio pulsars. 



3 EVIDENCE FOR LOW-MASS MILLISECOND PULSARS 

There are 1394 observational pulsars^ with known P and P (thus E) simultaneously. The 

numbers of millisecond and normal radio pulsars are 87 and 1307 if the dividing lines is by an 

"apparent" magnetic field of Bq = 6.4 x 10^^ \/ PP G= 5 x lO^'' G as shown in left upper panel 
of Fig.l. Here we only consider 87 millisecond pulsars {B < Bq G). 
It is known that the numerical value of rotation energy loss rate 

i; = -/nn:.3.ixio-(^f^p-3, (lo) 



^ http://www.atnf.csiro.au/research/pulsar/psrcat/ 



4 



Cui, Yue, Xu & Qiao 



Table 1 The radius and mass ratios of Group I to Group II from the relation of 
P — E and that of Group F to Group IF from the relation oi P — P. 



Relation 


Radius Radio 


Mass Ratio 


P- E 


0.32 


0.03 


P-P 


0.01 


0.13 



where rotational inertia of star / ~ |Afrg (the star is assumed to be a homogeneous rigid 
sphere), M = Mq, P_2o = P/10~^°. Note that the pulsar radius could be a variable. A 
potential drop across the open field line region for a rotating magnetic pulsar could be simply 

$ = 6.6 X 10^^Bi2rlp-^ V. (11) 

For giving the hints about the star radius and showing possible evidence for low-mass 
millisecond radio pulsars in the bare strange star model by observational data, we investigate 
the relations oi E — P and P — P as shown in left lower panel and right upper panel of Fig.l. 
Secondly, applying the Kmeans method, we divide them into two groups: 60 pulsars in Group I 

and 27 pulsars in Group 11; 63 in Group /'and 24 in Group IF. Thirdly, from B = 6.4x lO^^V PP 
G and assign P-20 = 1-0, $ = 1.0 x lO^^F in Eas. fTIHl and lfTT|l for 87 millisecond pulsars, we 
give the minimums of radius for these two groups rg/ ~ 0.35 and tqh sa 1.1 and the limit 
ranges of radii 0.065 < rg/' < 0.35 and 0.17 < re//' < 0.65. Finally, from the analysis above 
and according to Eq.(lO) in Xu (2005), we can obtain the radius and mass ratios of Group I 
to Group II from the relation oi P — E and those of Group I' to Group IF from the relation of 
P — P as shown in Table 1. 

4 CONCLUSIONS 

The total charge of an uniformly magnetized aligned pulsar with Goldreich- Julian density out- 
side the star is calculated. The boundary which separates the positive and negative flows in a 
magnetosphere is not necessary the critical lines if the real charge is not that given by Xu et 
al. (2006). The statistical relation between the pulsar's E (or the period derivative P) and the 
period P may indicate that the millisecond radio pulsars with small periods could be low-mass 
bare strange stars. 
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